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Abstract 

The measurement of — > t^v t at the B factories provides important constraints on the 
parameter tan (3/m H ± in the context of models with two Higgs doublets. Limits on this decay 
from e + e~ collisions at the Z peak were sensitive to the sum of B^ — ► r^v T and B^r — ► r^v T . Due 
to the possibly sizeable contribution from B^r — > t^v t we suggest that a signal for this combination 
might be observed if the LEP L3 Collaboration used their total data of ~ 3.6 x 10 6 hadronic decays 
of the Z boson. Moreover, we point out that a future Linear Collider operating at the Z peak (Giga 
Z option) could constrain tan/3/mjj± from the sum of these processes with a precision comparable 
to that anticipated at proposed high luminosity B factories from B — ► t^v t alone. 

PACS numbers: 12.60Fr, 13.20Hc 



I. INTRODUCTION 



In April 2006 the BELLE collaboration announced the first observation of the purely 
leptonic decay JB ± — > t^v t [l[ utilizing an integrated luminosity of 414 fb _1 . The measured 
branching ratio (BR) is in agreement with the Standard Model (SM) rate within theoretical 
and experimental errors: 

BR(5 ± - tV) = (1-79 ±° H(stat) t° t(syst)) x 10" 4 (1) 

Subsequently, the BABAR collaboration reported a measurement with an integrated lu- 
minosity of 346 fb -1 which is an average of separate analyses with semi-leptonic (2[ and 
hadronic jl] tags: 

BR(B ± -> r ± ^ T ) = (1.2 ± 0.4 ± 0.3 ± 0.2) x 10" 4 (2) 

The average of the BELLE and BABAR measurements is [3]: 

BR(5 ± -> t±v t ) = (1.41 t° .tl) x 10~ 4 (3) 

Significantly improved precision for BR( J B ± — > r 1 * 1 ^) would require a high luminosity C > 
10 35 cm -2 s -i B factory g g g g §, [Io|, [ll|. In the context of the SM the decay B ± -> r ± i/ T 
provides a direct measurement of the combination fsVub, where Jb is the decay constant 
which can only be calculated by non-perturbative techniques such as lattice QCD. Charged 
Higgs bosons [H^] present in the Two Higgs Doublet Model (2HDM) and the Minimal 
Supersymmetric SM (MSSM) would also mediate B ± — > t^v t [l2j with the New Physics 
contribution being sizeably enhanced if tan j3 (the ratio of vacuum expectation values of the 
two Higgs doublets) is large [l3j]. The above measurements of — > r ± z/ r now provide a 
very important constraint on the parameter tan/3/m#± in the context of the 2HDM and 
the MSSM. Hence this decay is of much interest in both the SM and models beyond the SM 
and improved precision in the above measurements is certainly desirable. 

Prior to the era of the B factories three LEP collaborations searched for B ± — > t^Vt 
and obtained upper bounds within an order of magnitude of the SM prediction [l4l 15. 16 



Such limits were actually sensitive to the sum of B ± — > t^v t and Bf — > r^v T [17J since 
the centre-of-mass energy (a/s = 91 GeV) was above the Bf production threshold (unlike 
the B factories). The strongest limits were set by the L3 collaboration which obtained 



BR(5 ± — > T^Ur) < 5.7 x 10~ 4 Uj. Since BR(5 ± — > t^Vt) has now been measured at the 
B factories, the L3 limit can now be used to provide a limit on the product of the transition 
probability f(b —> B c ) and BR(£>^ — > t^Vt). A quantitative study of the magnitude of 
the contribution of Bf — > t^v t to the LEP limits was performed in [13]. We update 
this analysis using the significant improvements in the measurements of the CKM matrix 
and calculations of Moreover, the measurements of the Bf production cross-section 
at the Fermilab Tevatron [HI, 0, 0, ED provide the first measurements of the transition 



probability for b — > B c and suggest much larger values than the theoretical estimations used 
in the numerical analysis of [17j . The L3 limit on BR(i?^ — > r^v^) was obtained with 
~ 1.5 x 10 6 hadronic Z boson decays, which is slightly less than half the full L3 data taken 
at the Z peak. We suggest that a search for B ± / Bf — > t^Vt- using the full L3 data sample 
(~ 3.6 x 10 6 hadronic Z decays [22|]) would not only strengthen the limit on the product 
of f{b — > B c ) and BR(i?^ — > t^Vt) but also offer the possibility of a signal, which would 
be an additional observation of B ± — >• t^v t and the first observation of B„ — > r^v T . It 
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G F = 1.16639 • 10" 5 GeV- 2 
= U.lUbDO Lrev 


m e = 0.511 MeV 

711 T = 1. ( ( ( Kje\ 


\V ub \ = 0.00386(28) 


\Vcb\ = 0.0416(9) 


m Bu = 5.279 GeV 
m Bc = 6.271 GeV 
f Bu = 0.216(22) GeV 


TBu = 1.638 • 10~ 12 s 
r Br = 0.463(177) • 10~ 12 s 
f Bc = 0.450 GeV 



TABLE I: Input parameters used in this paper, unless indicated otherwise in the text. 





BR(B+ -> t + v t ) 


BR(B+ 




BR(B+ e+v e ) 




1.2 • 10" 4 


5.5 


10- 7 


1.3 • 10~ n 


B c 


0.022 


9.3 


10~ 5 


2.2 • 10~ 9 



TABLE II: Standard Model predictions for the branching ratios (central values). 



is also pointed out that a future e+e Linear Collider operating at the Z peak (the Giga 

j) could offer similar sensitivity to the parameter tan/3/m#± from 



Z option [23j, |2J, [25 



these leptonic decays as the proposed high luminosity B factories. This article is structured 
as follows: in section II we present basic formulae for the decay rates for 5+/.B+ 



and discuss the H contribution; we study the admixture of Bf 
the Z peak in section III and give our conclusions in section IV. 



t+z/ and B ± — >• r+f at 



II. THE DECAYS B ± -> t±v AND B± -> t+1^ 

In the SM, the purely leptonic decays (^+i^) of B and -B+ proceed via annihilation to 
a W boson in the s-channel. The decay rate is given by (where g = uorc): 

r(fi, - i 1*) - iv^l [i-^rj W 

Due to helicity suppression, the rate is proportional to mj and one expects: 

BR(B+ -> r+z/ r ) : 55(5+ -> /x+^) : 5/2(5+ -> e+^ e ) = m 2 : mj : m 2 (5) 

These decays are relatively much more important for 5+ than 5+ due to the enhancement 
factor \V c b/V u b\ 2 (fB c / fB u ) 2 - Using the input parameters given in table [I], we obtain the SM 
predictions listed in table HT1 

The effect of 5+ in the 2HDM (Model II) on the decays 5+ -> £ + is e was considered in 
13| and the analogous analysis for 5+ —>■ l + vg was presented in 27]. In both cases the H ± 



contribution modifies the SM prediction by a global factor r q H where: 

r% = [l- tan 2 /? {M Bq /m H± ff = [1 - 5 2 M 2 J 2 (6) 

The 5 contribution interferes destructively with that of W . There are two solutions for 
r q H = 1 which occur at R = and R ~ 0.27 GeV" 1 for 5+ -> £+z/ £ (5 = and R ~ 0.26 
GeV -1 for 5+ — > This is shown in Fig. [1] where BR(B U — > t + u t ) is plotted as a 
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FIG. 1: BR{B U — > t + v t ) as a function of tan (3 and m^± . The plotted range of the BR corresponds 
to the l-o" range of the world average measurement (1.42 ± 0.43) • 10 -4 , and the line indicates the 
central value. 



function of tan/3 and m#±. For tan/3/m^i ~ 0, the BR remains at its SM value (slightly 
higher than the thin line indicating the central value of the experimental measurement), but 
this SM value can also be achieved along a line through the steep part of the surface where 
r„ = (l-2) 2 = l. 

If the b quark couples to both Higgs doublets at tree-level (which is referred to as the 



type-Ill 2HDM), Eq. © is modified to [28 



1 - 



tan 2 (3 



m 



B 



1 + e tan (3 m H± 



(7) 



In the Minimal Supersymmetric SM (MSSM), the parameter eo does not appear at tree- level 
but is generated at the 1-loop level (291 . 30 ] (with the main contribution originating from 
gluino diagrams) and may reach values of 0.01. The redefinition of both the b quark Yukawa 
coupling and the CKM matrix element V„h are encoded in i j3ll . 32| . The impact of e 7^ 
on rn has been developed in 33, 34, 35, 36^ . In particular, the value of R where r# = 1 
shifts depending on the magnitude and sign of e . 

In Fig. [2] we show the impact of the measurement of BR(B U — > t + u t ) on the plane of 
[tan j3, m H ±] in the 2HDM (Type III) which u pda tes the study of [28| (for a recent analogous 
plot with a somewhat lower value of fs see |37|). The white regions are excluded and the 
shaded areas correspond to BR{B U — > t + u t ) within the l-a experimental range. We plot 
overlapping bands for the l-a ranges of the input parameters and consider eo = 0, 0.01, —0.01. 
In the MSSM, positive values for e , corresponding to positive values of the /x-parameter, 
are preferred in order to explain the (g — 2) M anomaly [381 ] . but in general both signs are 
possible. 
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The different values for eo result in significantly different allowed regions in the plane 
of [tan/3,m#±]. Importantly, these constraints from BR(B U — > t + u t ) are from a tree- 
level process and when applied to the MSSM are only sensitive to the assumptions for the 
soft SUSY breaking sector via e (recently emphasized in (39J), i.e., a higher order effect. 
In contrast, other important B physics observables such as b — > 57, B s — B s mixing and 
Bd,s —> f^fi are & 11 1°°P induced processes. Consequently, constraints on the plane [tan (3,m H ±\ 
from such processes are very sensitive to the assumptions made for the sparticle masses, and 
in certain cases the constraints can be removed completely 1 . In global studies of B physics 

42| the measured BR(B U 



observables in specific MSSM scenarios [34J, |42[ the measured BR(B U — > r + z/ T ) also plays 
an important role. Certainly, improved precision for BR(B U — > t + v t ) is desirable and very 
relevant in the era of the LHC in which the plane [tan/3,m#±] will be probed via direct 
production of Higgs bosons. Currently only High Luminosity B factories operating at the 
T(4<S r ) are discussed when considering future facilities which could offer improved precision 
for BR(B U -> r + is T ). 

Another promising approach to probe the plan e \t an P ; rriff±] is via the tree-level 



MM, mm 



We note here that 
H, H and Df 
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contribution to the semileptonic decays B — > Drv 
can mediate the analogous leptonic decays K ± 
but constraints on the plane [tan/3,m#±] from these processes are not yet competitive. 
However, such processes might play a role in the future with increased experimental precision 
and reduced theoretical uncertainties. 



III. AT THE Z PEAK 



In this section we discuss the searches for B ± — > r^v using data from e + e~ collisions at 
the Z peak (y/s ~ 91 GeV). It was pointed out in 17j that such searches would also be 
sensitive to the decay B Assuming that the detection efficiencies are the same 

the ratio of r^v events originating from B ± — >• t^v and Bf — > r ± z/ is given by: 



Nc 



V, 



cb 



V, 



ub 



fib - S±) 



fB 



M Bc r Bc 



Mi 



(8) 



The largest uncertainty in the determination of N c is from the transition probability /(& — > 
Bf) and the decay constant fs c - The magnitude of N c is suppressed by the small f(b — ► Bf) 
but this can be compensated by the large ratio (V r hf r„) 2 / iVuhf r) 2 ■ Consequently N c can 
be similar in magnitude to N u . In the analysis of [17| three scenarios were defined in order 
to account for the error in the determination of N c /N u : "Central" and "Max/Min" (±lcr 
above/below the central values of the input parameters). Since the analysis of [l7[ there 
have been significant improvements in the measurements of V u b and V c b- In addition, the 
decay constant fs has now been calculated in unquenched lattice QCD with smaller errors 



and a central value considerably larger [46j] than the values used in both [17] and the L3 
analysis [3]. We are unaware of an unquenched lattice QCD calculation of fs c and the 



1 In the non-SUSY 2HDM (Model II) b — > sj constrains m H ± independently of tan j3. A recent study 
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4 if obtains m H ± > 295 GeV. 



In practice, the shorter lifetime of Bf would result in a slightly inferior detection efficiency 17 1 
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300 400 600 800 J 000 

FIG. 2: The constraint on the tan/3-m#± plane in the 2HDM (Type III) from the measurement of 
BR(B^ — > t^v t ). The coloured areas correspond to allowed ranges of tan/? and m H ± for various 
values of eo = 0,0.01,-0.01 (green, blue and red, respectively), BR(B U — ► t + i> t ) (l-er range, 
overlapping) and /b (overlapping). 



error in this parameter has not been reduced significantly since 17|. The main uncertainty 
in the ratio N c /N u is from f(b—> Bf), which in [13] was varied in the range suggested by 
theoretical estimations 0j: 2 x 10" 4 < f(b -> Bf) < 1 x 10~ 3 . At that time Bf was still 
undiscovered and hence there was no measurement of f(b — > Bf). 

However, f(b — > Bf) can now be extracted (although with a large uncertainty) from the 
measurement of ratio of Bf — > Jj^t^vg to B^ — > J/^>K ± which is defined by: 

_ a{Bj) ■ BR(B C -» J/^u e ) 
£ a(B+) ■ BR(B -> J/4)K+) 1 ' 



Tevatron Run II data gives lZ e = 0.28 ± 0.07 [20j, and the denominator in eq. Q) has been 
measured precisely by various experiments. The transition probability f(b — > B c ) determines 
cr(B^) and several theoretical calculations are available for BR(5 C — > J/^t^vg)). In Fig. [3] 
we display contours of lZ e as a function of BR(£> C — > J/\l/e + z/ e ) and f(b — > B c ), and the 
band denotes the prediction of the various theoretical calculations for BR(i? c — > J/^/e + u e ) 
whose values lie in the range (2.0 ~ 2.5)%_|48j. From Fig. [3] one can see that the Tevatron 
Run I measurement of lZ e = 0.13 ± 0.05 [l8[] is accommodated by f(b — > B c ) = 1.3 x 10 -3 . 
However, in order to satisfy the central value of the Run II measurement, the transition 
probability f(b —>■ B c ) needs to be 4.5 x 10 -3 . An even larger value for f(b —> B c ) was 
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4 5 6 

FIG. 3: Contours of lZ e in the plane of BR(B C — ► J/^e + u e ) and transition probability f(b — ► B c ). 
The shaded region denotes the theoretical prediction for BR(i? c — > J/^e^u e ). 



suggested in Ref. 49[. Such unexpectedly large values of /(6 — > £> c ), which are indicated 
by Tevatron Run II data, would significantly enhance the contribution of — > r ± i/ to the 
LEP searches for S 1 * 1 — > r^v. Of course, f(b — > 5 C ) is dependent on the available centre of 
mass energy (at higher energies there is more phase space to produce a charm quark instead 
of a light quark), but the value of f(b — > B c ) is expected to be of comparable size at LEP 
and at the LHC j47|. In our numerical analysis in Sec. IIII A\ we will consider values of 
f(b -> Bf) up to 5 ■ 10" 3 . 



A. The LEP search for B^ 1 — ► t^u and the contribution of B^r 



Three LEP collaborations searched for the decay JB ± — > t^v using data taken at the Z 
peak (y/s = 91 GeV). L3 14j used around 1.5 x 10 6 hadronic decays of the Z boson which 



corresponds to about half their total data [22|]. DELPHI 15] and ALEPH 16] used their 
full data samples of around 3.6 x 10 6 hadronic decays of the Z boson. The best sensitivity 
was from the L3 experiment which set the upper limit BR(5 ± — > r ± z/) < 5.7 x 10 -4 . The L3 
limit is of particular interest since it could be improved if the full data sample of ~ 3.6 x 10 6 
hadronic Z boson decays were used. 

The LEP searches were sensitive to events originating from both B and 
Bf — > t^v. Hence the published limits constrain the "effective branching ratio" defined by: 



BR 



eff 



BR(5 = 



N c 

1 + — 



(10) 



t^v at the Z peak. For searches at the T(4S) 
In our numerical analysis in this section we will 



This expression applies to searches for B 
clearly N c = and BR eff =BR(5 ± -> r ± z/) 
use eq. (11 01) with the experimental value of BR(S ± — > r^v) as input. The calculation of 
N c /N u in eq. (fTU]) uses eq .© (i.e. the expression for the SM) with input parameters taken 
from Table [H Our analysis can be applied to any model for which N c /N u ~ \N C / N u \sm, 
which includes the 2HDM because the scale factors in eq. ([6]) are almost equal. 
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The "max" scenario of [17j showed that the current limit of BR c fj < 5.7 x 10~ 4 would be 
sensitive to the SM rate for BK(B ± — * r^v). The measurements of BR(5 ± — > r ± z/) at the 
B factories are consistent with the SM prediction which suggests that the L3 search was not 
so far from observing a signal. 

In Fig. H] we plot BR off in the plane [f{b — > Bf), f Bc ], for two values of BR(5 ± — > r^v) 
corresponding to the central value and la above the world average. All other parameters 
are held at their central values from Tab. [B The region above the contour BR e gf = 5.7 x 10 -4 
(red/dark grey) is excluded by the L3 limit while the contour BR fj = 4x 10~ 4 represents 
the hypothetical sensitivity if the full data of 3.6 x 10 6 hadronic decays of the Z boson 
were used. The green/light grey area between the two contours is the area where a signal 
would be seen if the full dataset were studied. Depending on the other input parameters 
and the B — > r^v branching ratio, this area can cover a very significant part of the 
[f(b — > B^),f Bc ] parameter space. We therefore consider a re-analysis using the full L3 
dataset very worthwhile. 

A different way of studying the number of B c events was followed in [l7j]. The number 
of B c events per B u event can be calculated as a function of f(b — > Bf), and the authors 
obtained N c /N u = 1.2/(6 — > B^)/10~ 3 for central values of the input parameters ("max" 
scenario: 2.3). With updated values for the input parameters, we now find 



Nc _ 0.48 -/U^/IO" 3 (central values) 
N u " 1.50 -/U^/IO" 3 ("optimistic" values) 



3 ru^^„ (11) 



where for the "optimistic" values of the parameters from Tab. [I] we have chosen that end of 
the 1-cr range that results in a higher value for N c /N u , and the "optimistic" f Bc was chosen 
to be 550 MeV. 



These numbers are lower than those of [17J mainly because the central value of V u b/V c b 



has increased in the last ten years. The inverse of this ratio enters N c /N u quadratically and 
therefore reduces this quantity. On the other hand, experimental data does not preclude 
values of f(b — > Bf) which are much higher (a few xlO -3 ) than the theoretical estimates, 
and so the admixture of Bf — > r ± z/ can still easily reach 100%. 



B. Giga Z option at a future e+e Linear Collider 

A future e + e~ Linear Collider operating at the Z peak with a luminosity of 5 x 
10 33 cm~ 2 s _1 could produce 10 9 Z bosons in 50 — 100 days of operation 23, 24, 25, [26| ) 



This corresponds to roughly 1000 times the number of Z bosons recorded at each LEP de- 
tector. Historically, limits on B^ — * r^v from Z decays have been comparable to (if not 
stronger than) those at T(4S') for the same number of Z bosons and B mesons. For example, 
the CLEO collaboration obtained BR(B ± -> r ± i/) < 8.4 x 10" 4 with 9.7 x 10 6 B mesons 
while L3 obtained BR(B ± — > ^v) < 5.7 x 10~ 4 with 1.5 x 10 6 hadronic decays of the 



Z boson. 

High luminosity B factories [1, @, 0, S, @, 0, 11 1 anticipate data samples of 10 10 B mesons. 



By the time of operation of a Giga Z the two main sources of uncertainty in iV c (and hence 
BR e ff) will have been substantially reduced. The error in f(b — > B c ) will be reduced from 
LHC-b measurements jUJ of the cross-section in Eq. [91 and improved lattice calculations of 
f Bc and/or (/s c //s) would also reduce the error in iV c . In Table HTH we present the required 
number of B mesons and Z bosons for a precision of 20% and 4% in the measurement of 
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FIG. 4: The effective BR(B ± -> t±v) at the Z peak in the plane [f(b -> Bf), f Bc \- The published 
L3 limit and a possible stricter limit are indicated. In the upper (lower) panel BR(S =I: — > t^v) is 
taken to be the central value of the world average measurement (the 1-a upper value). 
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Error BR(B ± /B^ -> t±v) 


High Lum. i? Factory [B mesons) 


Giga Z (Z bosons) 


20% 
4% 


2.2 x 10 9 
8.1 x 10 10 


3.2 x 10 7 
8 x 10 8 



TABLE III: Required number of B mesons (Z bosons) for a precision of 20% and 4% in the 
measurement of BR(i? ± /i?^ — > r^u), assuming a signal of BR c g = 4 ± 2 x 10 -4 at L3. 



— > t^v at a high luminosity i? factory and BR fj at Giga Z . The numbers for a high 
luminosity B factory are taken from [1CJ |. For the Giga Z precision we assume a signal of 
BRefj = 4±2 x 10~ 4 (50% error) at L3 with 3.6 x 10 6 hadronic Z decays, and scale the error 
by l/y/N, where N is the total number of Z bosons at Giga Z divided by the full L3 data 
sample of ~ 5.1 x 10 6 Z bosons. 

It is clear from Table II I II that a Giga Z facility might be capable of measuring BR c g- 
in Eq. [10] with similar precision to that anticipated for B ± — > t^v at high luminosity B 
factories. We believe that this competitiveness of the Giga Z facility has not been pointed 
out for the leptonic B decays although it has been emphasized for the decay B — > X s vv 



in 



24i | . If both facilities were realized this would enable competitive and complementary 
constraints on tan/3/m#± in the context of models with H ± . 



IV. CONCLUSIONS 



The decay B ± — > r^v has been observed at the e + e~ B factories and is recognized as an 
important constraint on the parameter tan/3/m#± in the context of models with Two Higgs 
doublets. We studied the contribution of B„ —* t^v to the LEP searches for B^ — » r^v 



(first pointed out in [17j), whose main uncertainty is from the value for the transition 
probability b — > Bf which is now being measured at the Tevatron Run II. Using values 
of this transition probability which are consistent with the current Tevatron measurements 
(which accommodate values significantly larger than the theoretical estimations), we found 
that the contribution of Bf mesons to the search for B^ / Bf — > t^v can be as large as that 
of B ± . We suggested that a re-analysis of the L3 search for B^ — > r ± zy [lij using all the data 
taken at the Z peak could provide a signal for the admixture of B^ / Bf — > t^v. Finally, it 
was pointed out that the Giga Z option of a future e + e~ collider could offer measurements 
of these leptonic B ± /Bf decays which are comparable in precision and complementary with 
those anticipated at the proposed high luminosity B factories. 
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